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Abstract. Energies of graphene nanocones with 1 to 5 pentagonal disclinations are 
studied on an atomically detailed level. The numerical results are interpreted in terms 
of three different contributions to the cone energy: the core disclination energy, the 
bending energy of the cone surface, and the "line tension" energy of the cone edge 
that is related to different coordination of carbon atoms situated at the edge. This 
continuum description allows for a construction of analytic expressions for the cone 
energetics and indicates different regimes of cone sizes in which cones with a particular 
number of disclinations are preferred energywise. An important result of the study 
is that the energetics of various types of cones profoundly depends upon whether the 
dangling carbon bonds at the cone basis are saturated by hydrogen atoms or not. 
This may be of use for explaining the differences in the yields of various cone types in 
different production processes. 
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1. Introduction 

Carbon (C) atoms are famous for their ability to form various hollow structures pj, 
the best known of them being fullerenes[2] and carbon nanotubes[3]. Soon after the 
clear experimental detection of carbon nanotubes [3], conical carbon structures have 
been observed in the process of quenching hot carbon vapour on the graphite substrate 
[H [5]. These cones all had the same opening angle, the smallest one of the five that 
can be constructed from the honeycomb graphene structure (see below). It was later 
demonstrated that the structures of the same type can also be obtained by other 
experimental procedures, such as pyrolysis of hydrocarbons in plasma torch [6j. In this 
procedure, five different types of cone geometries were observed. The yields of conical 
structures in the prepared samples can be quite high (~ 20 %) [6]. This fact alone 
indicates the importance of such structures. In addition, conical graphene structures 
are also predicted to have very specific mechanical [7] and electronic [8j properties, 
and could possibly be used as nanoscopic electron-field emitters. Interestingly enough, 
graphite cones were found in natural samples |9] which additionally highlights the need 
to understand the occurrence of conical shapes made of carbon. 

The polyformity that is characteristic of hollow carbon shapes has its parallel in 
the structures formed by viral proteins. These also form hollow shells (the so called 
capsids) that protect the viral genetic material (DNA or RNA). Icosahedral viruses 
have a structure that is the same as the one exhibited by icosahedral fuUerenes |10j . 
The same proteins that form icosahedral shells may, under different conditions, also 
form tubular structures, similar in structure to carbon nanotubes. There seems to be a 
parallel for conical structures also - the core capsids of HIV viruses very much resemble 
(closed) graphene cone structures which was emphasised in Ref. [H]. The investigation 
of graphene cones may thus have an added value due to the parallehsm of structures 
found in graphene shapes and viral capsids. 

Unlike the fuUerenes and carbon nanotubes, conical carbon structures are relatively 
poorly explored. Not much is known about the optimal conditions for their production. 
A complete understanding of growth of such structures is also lacking, which is not too 
surprising since the details of growth of the much more investigated fuUerenes and 
carbon nanotubes are still missing [12]. Although there have been some proposals 
that the formation of cones is most likely due to their " designability" , i.e. the fact 
that they are a combinatorially favorable outcome of the growth [13], it is certainly 
important to first have a good understanding of and reliable estimates for the cone 
energies. Various contributions to the cone energetics are not completely self-evident, 
and it is not straightforward, for example, to estimate the importance of the curvature 
of the cone surface and its contribution to the total energy of the cone. 

In this article I shall examine the energetics of (single-walled) graphene cones. 
These cones contain multiple pentagonal defects. Cones with defective carbon rings 
whose number of sides is smaller than 5 are not considered in this study as these are 
shown to be less stable [l^. To understand various factors that influence the energy of 
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the cones numerical, all-atom calculations of the cone energetics are performed. The 
numerical results are interpreted in terms of different energy contributions that are 
disentangled from the total energy balance. This enabled a construction of reliable 
analytical expressions for the cone energetics. 

Section [2] discusses the details of the cone geometry and describes a particular 
construction scheme of the cones that is used in this study. 

Section [3] contains the details of the numerical calculations and their results. 

In Section H] the numerical data obtained in Sec. [3] is interpreted in a quite general 
framework that enables a construction of analytic expressions for the cone energies. 
Four contributions to the cone energy, each of them having a different physical origin, 
are identified. 

Section [5] contains a discussion of the results and some propositions that may relate 
the results of the article with the experimental data on carbon cones. 

2. Construction of cones 

There are five different classes of cones that can be constructed from the graphene lattice 
of carbon atoms. These differ in the total number of disclinations (i.e. pentagonal carbon 
rings) that they contain close to the cone apex, which is directly related to the cone 
opening angle, 9. The cones can be conceptually constructed by cutting out 60 degree 
wedges in the graphene plane, thereby creating a " cut-and-fold" cone "origami" that 
can be mathematically folded onto a three-dimensional pyramid-like object. Except 
for the cone with a single pentagonal disclination, this involves creasing along the 
lines connecting the neighboring pentagons. When the pyramid-like object is relaxed 
so that its energy is minimal, the shape acquires a cone-like appearance. Additional 
reconnections between the carbon atoms (e.g. between any two atoms that are on the 
cone edge and that could possibly form pentagons is they are sufficiently close [15]) are 
not allowed during the relaxation process. 

Except for the cone containing a single pentagonal ring close to its apex, all other 
cones can be constructed in a multitude of ways, depending on the arrangement of the 
pentagonal disclinations (see e.g. Ref. [IS]). These details are irrelevant for this study, 
as its main aim is to rationalize the cone energetics in the terms that belong to the 
realm of continuum physics. In each of the cone classes considered in the following, a 
fixed arrangement of the carbon pentagons is assumed. An isolated pentagon rule [17] 
was respected by all of the constructions so that the carbon pentagons do not share C-C 
bonds - each of the pentagons is completely surrounded by hexagons. The parameter 
that characterizes the cone is its size, i.e. the total number of carbon atoms that it 
contains. I shall also assume that the distances between the cone apex and its base, 
measured along the cone surface, are the same irrespective of the point on the base 
chosen. In other words, this means that the cones considered can be positioned on 
their base so that the cone apex projects exactly at the center of the base - the cones 
considered are right cones. 
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The five different cut-and-fold pieces of grapliene plane are illustrated in Fig. [TJ 
Similarly looking patterns were recently used for the construction of carbon nanotube 
caps, which necessarily have six pentagonal disclinations [H]. Note that the definition 
of the cone apex may be somewhat vague. In fact, it may be better to think of these 
structures as truncated cones or conical frusta. The smaller/top base of the conical 
frustum is then defined as the region surrounded the carbon pentagons (for graphene 
cone with a single disclination, the top base is the carbon pentagon itself). In the limit 
of very large cones, the distinction between the truncated cone and the regular cone 
becomes irrelevant since the radius of the bottom base of the cone is much larger than 
the radius of the cone top base. 

It is easy to show that the opening angle of the idealized cone [6, see Fig. [T]F)] can 
be related to the total number of disclinations (n) as 



This produces a discrete sequence of possible opening angles, 6'„ = 112.89, 83.62, 60, 
38.94, and 19.19 degrees, for n = 1, 5. It shall prove to be of use to consider a (planar) 
graphene disk as a special case of this construction. In this case, n = 0, and 9q = 180 
degrees. 

3. Numerical results for the cone energies 

The energetics of carbon-carbon bonding is simulated by using a second-generation 
reactive empirical bond-order potential by Brenner et al [19] which belongs to the Abell- 
Tersoff class of bond-order potentials. This potential is known to properly account for 
the anisotropy of C-C bonding. It also approximately includes the many-body effects. 
After the mathematical folding of the piece of graphene plane was performed and the 
broken C-C bonds reestablished along the edges of the shape (this part of the procedure 
also forms the carbon pentagons), the resulting structure is used as an initial guess for 
the conjugate gradient procedure that optimizes its shape so as to minimize the total 
energy. The numerical procedure used for this purpose is described in details in Ref. 
[20] . Similar computational procedure (the construction of initial mathematically folded 
shapes and subsequent relaxation using the conjugate gradient method in combination 
with Brenner's potential) has also been used in Ref. [2T] . 

Figure [2] displays the top and side views of the five types of cones that have 
been relaxed to their minimum energy configuration. The cones chosen for display 
are relatively small so that the details of the shape, including the pentagons around the 
cone apex, can be easily discerned. 

Figure [3] displays the energies of the cones depending on the number of C atoms 
they contain. In these simulations, the radius of the cut-out shape [denoted by R in 
Fig. [!][)] was varied. The energies displayed are in fact the excess energies (AE), i.e. 
the extra energy that the cones have with respect to the energy that the same number 
of carbon atoms would have in the infinite graphene plane. Brenner's potential [19] 




(1) 
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Figure 1. Construction of cones by creation of disclinations (triangular wedges/cuts) 
in the graphene lattice. The part of the graphene plane bounded by the thick lines is 
folded into a cone. Panels a), b), c), d), and e) display the graphene cut-outs used to 
construct the cones with one, two, three, four, and five carbon pentagons, respectively. 
The size of the cones is determined by the radius (i?) of the circle denoted by thin 
dashed line in panel e). Panel f) represents the cone and the quantities that define it 
(see text). 

predicts that the binding energy per carbon atom is = -7.39494 eV [2T], so that the 
energy of / = A^e^ was subtracted from the total energy of the cones in order to obtain 
the excess energy. The results for excess energies of graphene disks are also included in 
Fig. [3l Note the scatter of the data. It is not due to some numerical instability as shall 
be explained in the next section. 




Figure 2. Five different types of cones viewed from above (top row of images) and 
aside (bottom row of images). Panels a),b),c),d), and e) represent the n — 1, n = 2, 
n = 3, n = 4, and n = 5 cones, respectively. The cones contain 70, 70, 69, 72, and 74 
C atoms for n ~ 1, 2, 3, 4, and 5, respectively. 



4. General considerations of the cone energetics and interpretation of 
numerical data 

To understand the numerical results obtained in previous section, it helps to isolate four 
different contributions to the energy of a cone. First, there is a part of the energy that 
is directly related to the number of diclinations in the cone. Each of the pentagonal 
disclination (defect) carries the so-called disclination core energy (D) [211 122] that is 
related to change (reduction) of the local coordination of entities in the disclination. 
Second, the cones have a curved surface, so that there is a contribution to the total 
energy that is of the bending type (B) and directly related to nonvanishing curvature 
of the cone surface. For n this energy can be expressed as [2T] 

cos 6.fi ffmax\ 

B = TTCo . In , (2) 

Sm On V rmin J 

where Cq is the graphene bending rigidity [21], and r^m and r^ax are the distances 
from the cone apex to the approximate top and bottom bases of the cone, respectively 
[see Fig. [T]f)]. Equation [2] does not include the contribution related to bending of the 
shape surface in the top base, i.e. in the region where the disclinations accumulate. 
Assuming that the top base is nearly flat, or that the number of atoms that make it is 
much smaller from the total number of atoms in a whole conical shape, this part of the 
energy can be neglected. Third, there is an energy resulting from the coordination of 
the entities in a shape. This energy can be approximated as a sum of the energy that 
the same number of carbon atoms would have in the inflnite graphene sheet (negative, 
binding contribution, /) and the cone-edge energy related to reduced coordination of 
atoms located along the line that deflnes the bottom cone base (positive contribution, 
L). The effect of the reduced number of neighbors for the atom at the cone edge could 
also be called the edge-tension energy and is obviously proportional to the length of the 
cone base. The total cone energy, E is thus E = nD + B + I + L. Note that 5, /, and 
L depend on the size of the cone, i.e. on the total number of atoms, A^, that make it, 
but in a very different functional manner. The excess energy of the cones is obviously 
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Figure 3. Energies of carbon cones as a function of the total number of carbon 
atoms. In addition to results for n — 1, n — 2, n — 3, n — A, and n = 5 cones, the data 
for "n = cones" i.e. graphene disks is also shown. The numerical results in these 
plots are denoted by full squares. Full thick lines are a result of analytical expression 
for the cone energy, Eq. The thin dashed lines denote the bending (curvature) 

contribution to the cone energy, Eq. 



given by 

AE = uD + B + L. 
The total area, A of the conical shape can be approximately written as 



A 



. (en' 

TTSm - 



Or. 



sm 



(3) 



(4) 



so that the total number of atoms in the shape is 

N = A/ A,, (5) 

where the area per atom in the infinite planar sheet of graphene is given by Ac = rg-\/3/2. 
For the equilibrium nearest neighbor C-C distance in graphene, r^, Brenner's potential 
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predicts a value of [H] rg = 1.4204 A. The cone-edge energy depends on the total number 
of entities along the bottom basis of the cone. The total length of the bottom base edge 
is 2rmaxSin{9n/2)Tr, so that 

L = 26rmaxTrsm , (6) 

where S is the average line tension energy per unit length of the cone base. Combining 
equations ([3]), ([2]), and (jS]) yields the analytical expression for the dependence of the 
excess energy on the r^ax parameter of the cone: 

AE = nD + TTCo^^ In (^) + 25^. sin f ^) . (7) 

The dependence of r^ax on N can be obtained from equations (j4]) and ([5]) which yields 

(8) 



V r^ 



\ sin (^) 



mm 



sm 




As mentioned earlier, the r^m parameter is somewhat difficult to fix, but it can be 
estimated on the basis of constructions shown in Fig. [H Values of r^m = {ro/ sin(6'i/2), 
roy3/sin(^2/2), roVS/ sin(^3/2), roy21/[2 sin(^4/2)], 2roV3/ sin(^5/2)} were used for 
n = {1,2,3,4,5} constructions from Fig. [H respectively. For r^ax ^ '^mm, the precise 
value of rjjiin becomes irrelevant. It can be seen in Fig. [3] that Eq. ([71) (full thick 
lines) nicely accounts for numerical data on cone energetics for all n values. In these 
calculations, the core disclination energy and the graphene bending rigidity were set 
to D = 1.83 eV and Co=0.83 eV, respectively, in agreement with the values found in 
Ref. [21], and the best-fit value for S was found to be 5 = 1.30 eV/A (the same 
value was used for all n values of the cones). This should be the edge energy lost per 
unit length of the edge due to the reduced coordination of the carbon atoms situated 
on the cone edge. The atoms at the edge can have one, two, but also three nearest 
neighbors (see cones in [2]). The energy that is lost due to a missing neighbor should be 
about 2.46 eV (|e5|/3). The exact value of 6 depends on the edge of a particular cone 
considered, i.e. the number of singly, doubly, and triply coordinated carbon atoms and 
their distances. These characteristics depend both on the construction of the cone and 
on its total size (r^ax)- This explains the reason for scatter of the data in Fig. [3]- cones 
of different sizes generally have differently coordinated edges, so it may happen that in 
a cone of a particular size all of the edge atoms are doubly or triply coordinated, which 
is energetically favorable case (this is the case for all cones shown in Fig. [3]). Some 
other sizes may have edges that contain a lot of singly coordinated edge atoms. This 
is an interesting combinatorial problem that I shall not dwell too much upon, as the 
details of the shape that depend on the precise arrangement of atoms are not of interest 
to this article. Nevertheless, it should be noted that the value of 6 obtained from the 
fits can be nicely explained on the basis of energetics of C-C bonding. An important 
objection that could be made here is that the effective Brenner's potential cannot be 
reliably used in situations when there are dangling bonds. In Ref. [23] it was suggested 
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that Tersoff 's version of carbon-carbon bonding overestimates the excess danghng bond 
energy. Complex non-pairwise-additive effects may be also expected in cases when there 
are dangling bonds and these could be beyond the reach of the potential model used 
in this study. Nevertheless, these effects can be expected only to renormalize the value 
of 6, while the effective "line-tension" energy should still remain a useful concept even 
in that case. Even if the dangling bonds of C atoms at the cone edge were saturated 
by adsorption of hydrogen atoms [211 |25l [26], the reasoning in terms of the edge (line) 
tension energy would still be useful, only the value of 6 would differ (see below). It is 
also appropriate to note here that the calculations that treated the dangling bonds in 
the same fashion as described above have been used in a molecular dynamics scheme 
to simulate the dynamics of carbon "ribbons" that may form cylinders (nanotubes) by 
saturating some of their dangling bonds [T5] . 

Looking at the data in Fig. [3l it may seem that the edge tension energy dominates 
the cone energetics. This is evident both from the scatter of the data and from the 
apparent \/N behaviour of the excess energy. This is indeed the case except perhaps for 
n=5 cones, for which the bending contribution to the cone energy becomes important. 
This can be seen from the dashed lines in Fig. [3] that indicate the bending contribution 
to the cone energy [Eq. Q]. This is in fact one of the important results of this study: 
the bending/curvature contribution to the cone energetics is secondary when compared 
to the dangling bond/edge energy, so the considerations that concentrate exclusively 
on this part of the cone energy in order to explain the occurrence of the cones with 
particular opening angle in fact deal with rather small contribution to the total cone 
energy [1] (at least in cases when the bonds are not saturated by hydrogen atoms, see 
below). 

It is instructive to compare the analytic expressions for the cone energetics in a 
single plot, as is done in Fig. [H (note the log scales on both axes and that the curve 
corresponding to graphene disks is a straight line). It can be seen that the graphene 
disks have the smallest excess energies for small number of C atoms (A^ < 50). When 
the number of atoms in the structures reaches about = 50, the excess energies of all 
the cones except the n = 5 cone, become smaller than the excess energy of the disk. At 
this point the radius of the disk is about 0.7 nm. When A^ reaches about 150 atoms, 
the excess energy of n = 5 cones becomes the smallest and remains so as A^ increases. 
The transition region in the cone sizes that is observed in energetics of the cones smaller 
than about 150 atoms is a consequence of the interplay of different contributions to the 
excess energy of the cone. After a critical size of the graphene disk (A^ ^ 50) the cone 
structures become energetically preferable, since the energy that is lost due to creation 
of disclinations is compensated by the reduced length of the cone edges in comparison 
with the length of the graphene disk. For large enough cones (A^ > 150), this effect 
dominates the energetics of the structures considered, and n = 5 cones thus have the 
lowest excess energies, although they contain five pentagonal disclinations and their 
surface has the largest curvature. 
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Figure 4. Analytic (continuum) expressions for the cone energetics as a function of 
total number of atoms in the cones for the case when the dangling carbon bonds at 
the cone edge are not saturated (b = 1.30 eV/A). 

It should be instructive to see how all these conclusions change if the dangling bonds 
were saturated by hydrogen atoms. This may be the case for carbon cones obtained by 
pyrolysis of hydrocarbons [6]. The energy of C-C bond in graphene is[T9] 4.93 eV, while 
the energy of the C-H bond is [19] 4.526 eV. Thus the energy that is lost when replacing 
the C-C with C-H bond is 0.404 eV. This is also the energy that is lost per carbon 
atom at the cone edge when one of its bonds has been saturated by H atom. This is 
significantly smaller from 2.46 eV lost per carbon atom with a missing carbon neighbor. 
The line tension energy in this case is thus S = 0.404/2.46x1.30 eV/A= 0.213 eV/A, 
more than six times smaller than in the case of unsaturated bonds. Figure [5] shows 
how the energetics of cones and disks whose edges are saturated by hydrogen atoms 
changes with the total number of carbon atoms in the structures. It can be seen that 
the graphene disks are the structures with the lowest excess energies up to total number 
of about =2000 carbon atoms (compare this with =50 in the unsaturated bond 
case). When > 2000, the cones have lower excess energies than graphene disks. Again, 
the n = 5 cones are the most favourable in this respect. There is also an interesting 
crossover in cone energetics when A^ < 2000 - note how the curve for = 5 crosses 
the n = 4 and n = 3 curves for A^ ^ 170 and A^ ~ 700, respectively. In this case, the 
bending contribution to the cone energy becomes important even for large cones. For 
example, for n = 5 cone containing ?«9000 carbon atoms, the bending contribution to 
its energy is i? = 27.1 eV, while the edge energy is L = 39.1 eV. 




Figure 5. Analytic (continuum) expressions for the cone energetics as a function of 
total number of atoms in the cones for the case when the dangling carbon bonds at 
the cone edge are saturated with hydrogen atoms {6 — 0.213 eV/A). 



Intriguingly, for the case of carbon cones whose edge atoms are saturated by 
hydrogen atoms, the flat disk structures have the lowest excess energies up to a much 
larger number of atoms than in the case of unsaturated, dangling carbon bonds at the 
cone edges. The crossover region at which the cones become preferred energywise is in 
this case positioned at significantly larger number of carbon atoms (2000 vs 150 for the 
dangling bonds case). 

5. Discussion 

This study has demonstrated that graphene cones are superior to graphene disks with 
respect to their excess energy already for structures made of relatively small number of 
atoms (A^ > 50 for the dangling bond case and > 2000 for the hydrogen saturated 
bonds at the cone edge). For structures containing very large number of atoms, n = 5 
cones always have the smallest excess energy. For the dangling bond case, of all the 
shapes considered in this study that contain more than = 150 atoms, the n = 5 cones 
have the smallest excess energy. One could speculate that this may be one of the reasons 
for initial experimental observation of only n = 5 structures [HE]. Later experiments 
have detected all possible cones [6] and yields of cones of different n-type show a non- 
monotonic dependence on the number of disclinations, n. It would be tempting to assign 
this effect to the details of energetics of carbon cones in the cross-over region of sizes 
(50 < < 150). This explanation would require that the stable seeds for the cone 
growth are nucleated in this interesting region of sizes. Of course, one should keep in 
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mind that the growth of these structures is most hkely governed by kinetics and entropy 
in addition to energetics. 

Intriguingly, it has been shown that the energetics importantly changes when one 
considers the cones whose edges are saturated by hydrogen atoms. In that case, the 
crossover region is positioned at much larger number of atoms, and the flat graphene 
disks are the lowest excess energy structures up to about 2000 carbon atoms. It is 
thus plausible that the yields of different cone structures will be different depending on 
whether the hydrogen atoms are present in the production process, assuming of course 
that they have an active role in stabilizing the carbon edges. 

In summary, it has been demonstrated that the cone edge energy is the most 
important contribution to the total energy for cones that have more than about 150 
(2000) carbon atoms. The bending energy is significantly smaller from the dangling 
bond energy and becomes progressively less important for larger cones, since it scales 
with the total number of atoms, N, as log(iV), while the dangling bond energy scales as 
y/N. For cones whose edges are saturated by the hydrogen atoms, the consideration of 
the bending contribution to the total energy is necessary for the structures containing 
up to about 2000 carbon atoms. 
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